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Abstract 

Collective animal groups, the starling flocks particularly, may work at a 
critical state in which the efficiency of information propagation through the 
group is greatly increased. Here, we propose a minimum individual-based 
model to explore how the critical state can be reached. "Quorum response" , 
a type of social interaction in which an individual's chance of making one 
option is very likely if the number of its local neighbors committing to this 
option exceeds a threshold, is introduced as the local interaction. We high- 
light the potential for the enhanced ability of information transfer in the 
group at a critical state which endows each individual the ability to access 
information maximally to respond to environmental perturbations. The cor- 
relation function in 3D shows a scale free correlation which is supported by 
the observations in field experiments. 
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1. Introduction 

Collective animal groups, e.g. starling flocks, fish schools, are recently 
being called "collective minds" because of the coherent movements and 
the efficient responses to environmental perturbations as if the whole group 
is in one mind. It is obvious that the efficiency of information transfer in the 
group plays a key role in enabling these abilities. Cavagana et al. observed 
that in the airborne motion of large starling flocks, the length of correlation 
between two individuals' state scales with the size of the flock, the so called 
"scale-free correlation" [2J . This observation reveals that the starling flocks 
work at a critical state, in which one individual can effectively affect the 
state of any others' no matter what the group size is, and vice versa. This 
property confers the group an ability to share information efficiently so that 
it can optimally respond to external perturbations. A pioneering study on 
how information transfer in a collective animal group was carried out in a 

n 

fish school reacting to a risky perturbation in front [3]. It was found that the 
group members at the front made a quick rotation from the risk and their lo- 
cal neighbors behind imitated this behavior. The consecutive rotations of the 
group members resulted in a rapidly traveling "information waves", which 
rippled from the front to the rear at a speed 10 times faster than individual 
fish's speed. However, besides these experimental studies, the underlying 
micro-mechanism of information transfer is left largely ignored, particularly 
at an alarming situation in field [4J1. There are some other collective sys- 

n n 

terns, such as swarms of cancer cells 5 , bacterial colonies 6 and even human 

Q 

brains[7], share many similarities to the collective animal groups, and the 
working efficiency of which hugely depends on the underlying mechanism of 
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information transfer and processing. 

In this paper, we propose a minimum individual-based model to explore 
how the critical state of the collective animal groups, particularly the starling 
flocks in the airborne motion, can be reached. We apply an adapted form 
of "quorum response" as the local interaction rule and show that the group 
is poised at a critical state by tuning two parameters concerned with the 
distribution of individual's vigilance. The critical state is marked by a power 
law distribution of the size of the "information waves" -which is measured by 
counting the number of the informed individuals (turned from naive ones) 
in one run of information propagation. This feature endows each individual 
the potential ability to access information maximally to respond to environ- 
mental perturbations. We also calculate the correlation function in 3D and 
it shows a scale-free correlation which is supported by the observations in 
field experiments. 



2. Quorum response 



Quorum response is an interaction widely found in the bee and ant 



colonies 8, 9|, the cockroach aggregations 101, the broiler chicken crowds 111 

n 

and the fish schools [12J. The essence of quorum response is that an individ- 
ual is very likely to choose one option among others if the number of its local 
neighbors committing to this option exceeds a threshold. Let's consider the 
following simple example, suppose there are only two options, e.g. being at 
an alarmed state ("+" state) or a naive one ("-" state). The mathematical 
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description of quorum response is as follows 



p\(t + i)= , pi(t+i) = i- P %(t+i) (i) 



where p\(t + 1) is the probability of the individual i choosing to be at an 
alarmed state at time step t + 1, Uq is the constant total number of local 
neighbors to the individual i and n l + (t) or n l _(t) is the number of local neigh- 
bors who commits to the "+" or "-" option at time t, respectively. Obviously, 
the equation n\ (t)+ri i _(t) = no holds. The term n * «j plays the role of the 
threshold value for individual i. Parameters ctj, which quantifies how sensi- 
tive the individual i responds to its local neighbors' commitments is coined 
as "vigilance number" below. Parameter k determines the steepness of the 
function. In Figure 1, we see that the probability is a monotonic increasing 
function. Near the threshold value, it has an inflection point with a rapid 
increase and the function is sigmoid. When the parameter k approaches 10, 
the plot is practically a step-like switch at the threshold value, jumping from 
zero to unity. Quorum response is essentially a distributed positive feed- 
back process that enables information propagation. It is believed that this 
interaction can enhance decision speed and accuracy for a group to make a 
collective decision ljj. 

Until now, because of the complexity of the topic, how individuals in the 
collective animal group interact with each other in the field does not have a 



conclusive answer 
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151 ] . For simplicity, we adapt the interaction of quorum 
response as a step-like form function so that the focal individual's probability 
of taking an option is or 1 depending on wether the number of its nearest 
neighbors committing to this option exceeds the threshold value of the focal 
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Figure 1: Function of quorum response according to equation (1). The y axis is the 
probability for the individual i at time step t + 1 to choose the option (+) and the x 
axis is the number of the local neighbors who commit to the option (+) at time t. The 
total number of local neighbors to the individual i is set to 100, n = 100. Parameter 
k — 10 determines the steepness of the function, ai, which quantifies the sensitivity of the 
individual i to the commitments of its local neighbors, is coined as "vigilance number" to 
the individual i. It is related to the threshold value by n * on. 

individual. 

3. A minimum model 

A group is composed of iV individuals, in which each one is assigned 
a "vigilance number" ctj. We assume that the distribution of ctj obeys a 
truncated Gaussian distribution in the interval (0, 1), with the mean being 
a and standard deviation being A. Each individual is fixed on an evenly 
spaced grid during the process of information propagation. The threshold 
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value of the individual i is defined as, 

Ti = n * ati 

where no is the constant number of nearest neighbors to an individual (no = 
4 in 2D and no = 6 in 3D in the model, with some modifications to the 
individuals at boundaries). Each individual can either be in an alarmed 
state ("+" state) a = 1 or in a naive state ("-"state) a = — 1. The only rule 
of local interaction is the adapted quorum response rule: 

if Ti < n\(i) then a^t' > t) = 1 

If the condition is satisfied, then the individual i will turn into an alarmed 
one <Ti(t + 1) = 1 at the next time step and remains unchanged during the 
process of information propagation. If not satisfied, and if the individual i is 
at a naive state, then it will remain at a naive state <7j(i+l) = —1 at the next 
step. When the individual i turned from a naive state into an alarmed one, 
it will continue to interact with its nearest neighbors. The local interactions 
may eventually cascade into "information waves" at a global level. Initially 
all the individuals are set to naive state ( "-" state) and information is tapped 
into the system by randomly picking one individual at the peripheries and 
turning its state to the alarmed one ("+" state). A process of information 
transfer is considered completed when all the alarmed individuals don't have 
the ability to alter the state of its nearest neighbors anymore. 

We find that if the parameters of a and A are fine tuned- a = 0.28, 
A = 0.15 for 2D and a = 0.29, A = 0.35 for 3D to be specific- a power 
law distribution of the size of information waves is emerged (see Figure 2). 
The size of the information waves is quantified by counting the number of 
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Size of information waves 



Figure 2: Distribution of the size of the information wave, (a) 2D 100x100 array. The 
fine tuned parameters in a Gaussian distribution of on are mean a = 0.28 and standard 
deviation A = 0.15, and data is averaged over 5 x 10 5 simulation; (b) 3D 25x25x25 array. 
The parameters are mean a — 0.29 and standard deviation A = 0.35, and data is averaged 
over 5 x 10 5 simulation. Note the double-log of the plot. 

the affected individuals who turned from naive state into alarmed one in 
one run of information transfer process. The curve is linear for more than 
three decades, with slope of the fitting line is -1.27 for 2D and -1.31 for 3D, 
indicates that the fine-tuned parameters have poised the system to a critical 
state[ii|. Suppose the mean value a is lower, the whole group is poised 
at a super-critical state, where any small perturbations will very likely to 
cascade into big information waves that affect a lot of individuals. It will 
cause a lot of energy waste. If the mean value a is higher, the whole group is 
poised at a sub-critical state, the perturbations only affect in local areas, and 
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information transfer is greatly hampered. The collectively-tuned parameters 
of a and A, which poise the group to a critical state, are resulted from 
a long-term adaptation in the risky naturejl^]. In the model, for clarity, 
we assume the vigilance number of each individual in the group obeys a 
same Gaussian distribution. Yet it was observed in field that individuals 



at periphery are more vigi 
the so called "edge effects" 



ant in average than their conspecifics in center, 
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191 ]. This fact will increase the probability of 



the big information waves arising, which practically enhance communication 
between individuals j^ll]. 

We would like to make a comparison with the well-studied theory of self- 
organized criticality (SOC), which was proposed by Bak, Tang and Wiesen- 



201 ] . It is now a commonly accepted underlying mechanism to 



feld in 1987 

a wide range of phenomena such as earthquakes [23{], avalanches in granular 



materials 



16]etc, and even applied to the dynamics of neutral networks 25] 



and natural evolution 



241 ]. It states that a complex system can organize itself 



to a critical state in a dynamic process without any tuning of parameters 
from outside. The mark of a system entering a critical state is a power law 
distribution of the size of avalanches- the counterpart of information wave- 
which is very similar to our model. However, the theory of SOC states that 
the self-organized critical state, which is not necessarily subject to natural 
selection, is a robust attractor for a complex system no matter what the 
initial state is. On the contrary, the critical state of our model is already 
collectively-tuned under the selection pressure in the risky environments. 



8 



1.0 
0.8 
0.6 
0.4 
0.2 
0.0 



L=20 
L=22 
L=25 



0.0 0^2 0^4 0^6 0^8 l!o 

r/L 



Figure 3: Correlation Versus the rescaled length. Correlation function, at a critical state, 
of groups with different size is plotted versus the rescaled distance r/L. 

4. A scale-free correlation function 

There are already observations from field experiments of the correlation 
function of the individuals' states in starling flocks j^j]. The formula of the 
correlation functions is, 

= 1 Ejj(Tiaj5(r - r i3 ) 

which states how strength of the influence from one bird to any other one 
varies with the distance r. Delta function S(r— r^) picks a pair of individual i 
and j at a distance of r. cq is a normalization factor such that C(r — 0) = 1. 
In the field experiments of starling flocks, it was found there are two large 
clusters of strongly correlated birds in the flock, where individuals in each 
cluster have the same velocity state. In order to mimic this condition, we 
only calculate the correlation function when a sufficient large information 
wave (> 0.3 • N) is propagated so that the group is divided into two clusters 
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of alarmed individuals and naive ones at a global level. In each simulation 
50 individuals are randomly sampled and the correlation function is aver- 
aged among the sample. The distance was calculated in discrete integers, 

= sj (xi — Xj) 2 + (yi — ?/j) 2 J , where [J is the floor function. We find that 
when the system is tuned into a critical state the correlation function of the 
group at different size collapse together (see Figure 3), which shows that 
two individuals are correlated in a way that the correlation strength is not 
relevant to the group size, a scale-free correlation. 

Recently, many studies focus on how the ordered state is resulted from 
local interactions and how the directional information is propagated through 



the group. Theoretical 
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4| and experimental |27| studies have revealed 
that it is a universal property of a collective animal group to evolve into an 
ordered state by applying "averaging the velocity among its local neighbors" 
as the interaction rule. By this way the informed individual (s) has to interact 
with the naive one(s) repeatedly to spread the directional information j^]. It 
is a gradual process, which is very unlikely that the alarming information 
is propagated this way. On the contrary, applying the interaction rule of 
quorum response and positioning the system at a critical state, individuals 
can respond to the risky environments in a drastic way by cascading a local 
perturbation to a big information wave. Thus information is propagated 
more efficiently. 

5. Conclusions 

We proposed a minimum model, with introduction of an adapted form 
of quorum response as the local interaction rule among individuals, to study 



10 



how information is propagated in the collective animal group. By tuning the 
parameters of the mean and the standard deviation of a Gaussian distribution 
of individual's "vigilance number", we found that the group is poised to a 
critical state. It is reflected by a power law distribution of the size of the 
information waves. Being poised at a critical state, each individual can access 
information maximally to respond to environmental perturbations. We also 
calculated the correlation function of different sized groups in 3D and found 
that date collapsed on one curve. This demonstrates that the correlation is 
scale free, which is supported by the observations in field experiments. 

We hope this minimum model captures the essences of information trans- 
fer in the collective animal group and highlights the advantages of critical 
state underlying information transfer. Recently, it was found that the cortical 
neuron networks work at a critical state, on which the information transmis- 



sion, adaptability and computational power would be optimized 2l|, |25j. It 



seems the "collective minds" shares deep similarities to "minds" indeed. 
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